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'.THE COMPRESSIBLE POTENTIAL ELOW PAST ELLIPTIC 



SYMMETRICAL CYLINDERS AT ZERO ANGLE ' OE ATTACK 



AND WITH NO CIRCULATION* 



By W. Hantzsche and H. Wendt 



For the tunnel corrections of compressible flows 
those profile" are of interest for which at least the 
second approximation of the Janzen-Rayleigh method can be 
applied in clos.ed| form. One such case is presented by- 



certain ellips-oXd-al symmetrical cylinders located in the 
center of a tunnel with fixed walls and whose maximum ve- 
locity, incompressible, is twice the velocity of flow. 
In ..the numerical solution the maximum velocity at the 
profile and -the tunnel wall a's well as the' entry of sonic 
velo.city is computed. > - 

The. velocity distribution past, .the- - cont our and in 
the minimum cross section .at various Mach numbers is' il- 
lustrated on a worked out-example. 



The method of Janzen (reference l) and Rayleigh (ref- 
erence 2) for the step-by-step calculation of subsonic 
compressible potential flow from the incompressible flow 
was in original form applied only to flow past the circu- . 
lar profile and to the sphere in free air. Poggi modified 
the solution of the boundary value problem for two- 
dimensional flows so as to make the method equally appli- 
cable to other simple profiles (reference 3). The math- 
ematical task is largely confined to the determination of 
the velocity or pressure distribution on the profile. 
Eor the solution of the second step- (development of the 
potential to the quadratic term of the Mach number), Imai 

*"Die kompressible PotentialstrSmung urn eine Schar von 
nichtangestellten symmetri.schen Zylindern im. Kanal , " 
Luf tfahrtf o'rschung, vol. 18, no. 9., Sept,. 20, 1941, 
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and Aihari (reference 4) ha'ye indicated a method by intro- 
ducing the variables z an<| z instead of x and y 
previously employed "by Buseaann (reference 5) in his 
treatment of the incompressible flow, which enables the 
obtainment of the potential in the entire field even for 
more' complicated cases. By-way of application Imai and 
Aihari calculated the flow past an ellipse in free air. 
In tylSifiWtfsft 1 W*Xaplan ! s calculation patterned after Poggi 
(reference 6) for the ellipse at zero angle of attack and 
with no circulation, which requires the series, the solu- 
tion can he given in closed form. 

For flows in the tunnel only the second approximation 
about an almost circular contour by S. Lamia (reference 7) 
is available. He used the method of Jansen and Rayleigh, 
with the incompressible flow produced by superposition of 
a parallel flow with a doublet and its reflections at the 
tunnel walls as first approximation. 

The present report deals with the second step for 
the incompressible flow about a group of ellipsoidal cyl- 
inders in the tunnel. The profile curves are denoted by 
circles after conformal transformation of the strip 
bounded by two parallels, mapped on a plane secti onali zed 
along a semistraight line. In incompressible flow the 

cylinders have for each thickness ratio A hicknesg \ 

\ chord / 

double the flow velocity as maximum velocity. The dif- 
ferential equation is integrated according to Imai and 
Aihari . 



I. DIFFERENTIAL EQUATION OF THE FLOW . ' 

The potential $ (x, y) of a stationary, nonrota- 
tional two-dimensional compressible flow satisfies the 
differential equation of the second order:' 

where u = $ x , v =- $ y indicate , ' respectively , ' the com- 
ponents of the velocity' along the .x- and y-axis- of a 
rectangular system' of coordinates x,y and- a- the sonic 
velocity. This is a function of u,v; with a. m as sonic 

velocity and U as the flow velocity at infinity we get 



WACA Technical Memorandum No; 1030 



a 2 = : k 2 OT - (u 2 ' + v 2? - if T (2) 



Equation (l) can "be' written in' the form 

$ xx + *yy = ^kx^x 2 + - 2$ x $ y <5 xy + $ yy $y 2 ] / (3) 

Envisaging the potential, developed along the Mach number 
it = -2- 

3. CO . 

$ = $o'+ $! M 2 + . . . 
equation (3) affords differential equations of the form 

A<S 0 = 0 (4) 

A $1 = ^ ((a» 0 ) xx ($o)x 3+ 2^o) X ^o)y^o) X y + .^ 0 )yy^oV 2 } (5) 

$ 0 is the potential in incompressible flow; the deter- 
mination of $j represents the second step of the Janzen- 

Rayleigh method of approximation. As potential of an in- 
compressible flow $ 0 is the real part of an analytical 

, v f (z ) + f (z ) 
function of z = x + iy : 6 0 = Re f(z) = — '—— » if 

. . 3. 



f(z) is the -conjugate complex value of f(z). 

With r , 

, x _ ■ df 1 df ,. df 

' . • df _1_ f 'df df 1 

c$ 0 ; y = - Jm dz = 2i dz + dz j 

(<h ) = Re 



dz 




f 




2 




df. 


f 


dz 


a 


d a 


f 



/ N _ d 2 f i r d s f - df_fi 

^o'xy = " dz 3 = 2i l"'.d« a dzsj 
\ -p d 2 f 1 fd 2 f d 2 ?] 
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the right-hand side of equation (5) 'can be formed to read: 



i fd a f (MS\ S a. iff fiT\ 8 \ 

2TJ 3 \dz 3 Vdz / dz 3 Vdz/ J 
The equation for $ 1 assumes the form 

Introducing the complex variables z = x + iy and 

z = x - iy in place of x and y it can he transformed 

with and & i (a AN t0 

ox oz oz dy \oz dz/ 



dzoz " 4U 3 — \dz^ Uzy J 
the integral of which can be indicated 

•i - M {S / (if) 3 *•♦»<•> + » (7)} 

The analytical functions 3? (z) and G (z") are arbitrary, 
pending determination by the boundary conditions. 

The formal use of z and F as independent vari- 
ables still requires justification. 



The conjugate complex value f(z) of the analytical 
functio n f (z) is no longer an analytical function of z, 
while fTzT is obviously analytical in z". To express 
this we write 

■ f (z) = f(z) (8) 

Since dz = dz , the derivatives of f with respect to 
t become 
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df (z ) _ df (¥ ) , dgf (z ) _ d s f (z ) 
dz &T dz dz" 

whence equation (6) can "be written 



(9) 



«• ■ (ft)' ♦ B (^F)'] (l0) 

Permitting the independent variables x and y to 
assume temporarily any complex values Xj + ix 3 and 

yi + iy3f then f(z), f("z) and $ ostensibly "become 

analytical functions of the two complex variables x and 
y. Then the substitution of the independent variables 

Zi = x + iy 

and z s = x - iy 

for x and y changes equation (10) to 



* 3 $1 = 
d z ! d z 3 8U 



a 



d 3 f (z a 



dz 



.3 



> Lixi y + d3f iix) / df (z 3 ) \ s i 

\ dzj, y dzj 3 V dz 3 / J 



So for we get 



= _1^ I df(z ? ) /> /df(z a )N 3 
1 8U 3 L dz 3 J \ dz a / 1 

+ af(» 4 ) y ( d ^ Zg) ) a dz 3 + 23f + 2G(z 3 )J (11) 

with F (zj and 0 (z 3 ) indicating the analytical 

functions of Zj and z 3l ' respectively. Again limited 

to real x and y affords z x = z and z 3 = "ST. Taking 

as moreover the real part of equation (ll) so as to secure a 

real solution of equation ( 6 ) , the f inal expres si on 
(equation (7)) with due account of equations (8) and (9) 
reads : 

•i - & * &Pf (£)' d2 + r M * 8 m ] 

I I : 
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If the potential is referred to the complex coordi- 
nates £ and £ of a £-plane tied to the xy-plane by 
the conformal transformation t, = £(z), then $ x assumes 

the form 




Conformal transformation is usually resorted to for 
reasons of easier compliance with the boundary condition 
that the derivation of $ x normal to the profile contour 

should disappear. 

The differential equation for the second term of the 
stream function in the development along the Mach number 
can "be transformed in analogous manner, 

The stream function i|/(x,y) of a two-dimensional 
stationary nonrotat ional compressible flow meets the dif- 
ferential equation 

with up = p ro "4/y and vp = - \|/ x the components of 
the stream density. 

With \|f = \|f + M s + . . . (ty stream function of 
incompressible flow) equation (13) affords for \|/ a the 
differential equation 

^1 - <*o>xx - 2 ^o>xy U5 + ( ^o)yy . 

whence the introduction of z = x + iy and z" = x - iy , 
if \j/ Q = Jm_ f ( z ) , gives 

v 1 |T 4U 3 — \_ dz 3 \ dz / J 

or 



NACA Technical. Memorandum No.. .1030. 



7 



A comparison of equation (7) and equation (14) 
shows that the complex functions occurring at $ a and \|/.; 

behind the sign Re and Jm differ only in the terms 
added "by the boundary conditions . ■ 

II. INCOMPRESSIBLE FLOW IN THE TUNNEL 



We proceed from an arbitrary symmetrical profile at 
zero angle of attack placed in the center of a channel 
with fixed walls. The flow of the entire plane obtained 
by reflection at the tunnel walls is referred to a system 
of rectangular cartesian coordinates x, y, the x-axis 
of which is located at tunnel center. The tunnel height 
can be chosen equal to 2tt without limitation of gener- 
ality. If z = x + iy again indicates the complex co- 
ordinate in the xy-plane , the complex potential of the 
flow is an analytical function f(z), which at infinity 
acts as Uz (U = flow velocity) and otherwise has singu- 
larities only on the inside of the profile; f(z) is a 
simple periodic function with the period 2 tt i . Then the 
transformation £ 1 = e * maps the z-plane into the com- 
plex £' = I* + in'-piane and each strip parallel to the 
x-axis of width 2 tt covers the entire plane. Consider- 
ing the strip - tt ^ y ^ tt in particular, we find the 
tunnel walls y = - tt and y = tt changing to the twice 
traversed negative | r -axis. The line of symmetry of the 
tunnel becomes the positive |'-axis, the symmetrical 
profile changes to one symmetrical to the | x -axis. The 
infinitely remote point of the z-plane is reflected in 
the zero point and infinitely. remote .point of the £ '-plane 
Erom the behavior of . f(z)\ for z = oo it follows that in 
the vicinity of = 0 and £ T = t», f(£ l ) acts as... 

U In £. Aside from these two points the function has no 
singular- points outside of the profiles; f(£*) is there- 
fore the potential of the flow in the £' -plane", produced 
by a source with 'the' yield 2 tt U • in . the zero point. The 
profile and the ^ '-axis are streamlines. 

Concentrating on the group of profile contours in 
the tunnel which the transformation £ T = e z maps into 
circles, the complex potential f(£ ! ) is readily.. 
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indicated by ref lection . on' th.e. circle? 

f(f) = u m ^' ( ^f_V C) (15) 

if h is the center abs.cissa of the circle of radius a 
and c = t~ " 

D 

The family of profile curves in the t;unnel is shown 
in figure 1 a. 

These profile curve's are symmetrical with the median 
line of the tunnel and at right angle to it , as is readily 
seen from figure lb conformable to the chord-tangent 
theorem. The highest point of a profile curve in the 
£'-plane is given by the contact point of the tangent 
placed on the circle from the origin of the coordinates. 
Obviously the same profile curve is obtained for values 
of a and b with constant ratio j~ up to a parallel 

displacement along the tunnel axis. Thus keeping the 
center distance b fixed while varying the radius a 
from zero. to b-. affords the whole group; the position 
shown in figure 1 a then is' obtained by parallel shifting 
of the individual curves. The relationship of displace- 
ment ratio, to thickness ratio of the profile curves is' 
indicated in figure 2. "The maximum velocity follows 
from 

df _ df d£' _ df_ ? ,' 
dz • d£« dz ~, d£' 5 

for the contact- point of the. tangents 

i » = 7b 3 - a 3 ^ 133 ~ a * + i ^] 



*To find t'he flow represented by. f(£') if the profile 
in the £ '-plane is not yet a circle, the outer zone of 
the -profile would have to be mapped onto the. outer zone 
of a circle of a complex ■ z '-plane. Of the con-formal 
transformation it is demanded that it be symmetrical to 
the £ * -axis, that is, the £ 1 -axi s.. b ecome s ...the x'-axis 
and the infinitely remote point of the £'-plane is 
mapped in that of the z'-plane. Then the flow in the 
z'-plane again becomes the flow of a source about a circle. 
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with allowance for equation (15) at 

df . _ 

— = u - i v = 2U 
dz 

hence the maximum velocity is always twice the flow ve- 
locity, as is readily understood for the two extreme 
cases a—* 0 (sphere in free air) and a — » "b (very long 
plate of thickness ^ in tunnel of height h). 



III. COMPRESSIBLE FLOW IN THE TUNNEL 



In the following the circular plane is referred to 
a coordinate system £ = £ ' - b through the center of 
the circle. 

Then the complex potential (equation (15)) has the 

form 

fit) = u m ( £ + ^Mi + °> (is) 

whence the incompressible velocity potential 

$ 0 = M f(&) ■ He U In U * * °° . (17) 

The additional potential is according to equation 

(12) 1 

i ». - }- He f- M-±n_L + *± _ * in (1 + * > 
U *i 4 — L s u £ 

* ,„ & + ° + - c ) _ h(h ~,o) 

+ hd m -y— + ru t c) , £(Tt c) 

+ (b_- cjd ln . + (p - c2 In (t + h) 

■ I + c L - £ + c 

(h - c)_d ln (£ + c ) (b - c) 3 

£ + c " ""(j + o)Cj + o) 

+ Ij <£) + Gj. (J)) (18) 
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with d = - 1 - 2£± 

a 3 

Because of the conformal transformation of the 
z-plane into the £-plane the "boundary con- 

ditions -in the £-plane can he easily indicated. 

1. .. The disappearance of 'the normal component of the 
velocity on the profile of the xy-plane stipulates that 
the derivation of § 1 normal to the circle r = a in 

the £-plane he zero: 

' '(&)„-''•" 

2. The tunnel walls must he streamlines, hence the 

0$ 

derivation -r — on their corresponding part of the neg- 
ocp 

ative | -axis must disappear: 

r\ = o 



3. The additional velocity at infinity must "become 
zero in the plane of the tunnel. In ;point £ = - h and 
t = 00 the expressions 



and 



**i = ^ k it + h) + 2*L J, (J + o) 
ox 0 & dri 



- Jm (J + h) + Re (£ + h) 



must therefore disappear. 

4. As the flow is to "be without circulation, the 

•derivative for c = ±a must hecome zero: 

dcp * 
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The functions F a (£) 
and &! (£) ehould be so de 
fined that <!>! (£) complies 
with the conditions 1, 2, 3, 
and 4, and that the velocity 
field in the xy-plane has no 
singularities outside of the 
profile at finity 

It then affords for $j 



The procedure of estab- 
lishing 7 2 (£) and (f) 
is given elsewhere. The ex- / 0 
i stent integrals are i-nd-eter- "V 
ai*fftT integrals, since a con- L ' 
stant amounts to nothing at 
the potential. The integra- 
tions are easily carried out. 
They were omitted here "because 
it makes the expressions more 
complicated and is moreover 
unnecessary for the predic- 
tion of the velocity field. 



L C + c c(a a + c£) J (a 2 - 



' + (a 2 + c£) 

>ln(f + ft) 



df II 



.(l+8^_a^.)[ln(t + 4)-lnf + ln(t + «)]J Ill 



(19) 
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So far the "boundary condi- 
tions 1 hare "been net without 
taking the others into account. 
This was carried out separately 
for each term on the right-hand 
side of equation (18). 
on the right-hand side 
form Re. (h(£ ) k(f)) and 
nishes as constituent 

radial derivative -r- 



to 



fur- 
the 



^•*<C)^ + jfMC)e--] r = a 



This term is now 
part of a function of 

Re_ (g(e iqp )). Forming both 
pressions 



the individual terms of equa- 
tion (18). The two underscored 
terms in equation (18) together 
are purely imaginary, hence may 
"be omitted. These were in- 
One term eluded in equation (18) only 
has the because of the parallel to the 
stream function. 

The second boundary condi- 

^jf) - 0 i e satisfied 

T|«0 

by all nonlogarithmic terms, be- 
cause the expression below the 
sign He. derivated with re- 
spect to <p becomes purely im- 
aginary. Otherwise it may hap- 
pen, for instance, at term 



for 



tion 



(3 



a real 

e^, 



say, 
ex- 



- as 



51 



(9 



if 



and — 




(third line), that 



they patently have on circle 
r ■ a the radial derivative: 

- Re g (e icp ) precisely required 
for the compensation. 

The choice of either func- 
tion is, for the time being, 
open; a decision is usually ar- 
rived at because only one of the 
two functions satisfies the 
boundary conditions 2 at the 
same time. It therefore seemed 
advisable to treat the two addi- 
tive components of g (e*^) dif- 
ferently on many terms. 

By the present method of 
compliance with conditions 1 
both the real and imaginary parts 
of function g (e iq P) at the cir- 
cle are compensated, which is 
not at all necessary. 

The first eight lines of 
equation (19) indicate the pro- 
portion for the additive_f unc- 
tions 7a '<£) and Ox (£) for 



In U + b) 
t 

real proportions remain. To 
bring these to disappearance 
the two lines of equation (19) 
below boundary conditions 2 
were necessary. The expres- 
sions are so chosen that condi- 
tions 1 are always complied 
with. 

Jiastly, the solution of 
the homogeneous equation in the 
last line has been given a sui- 
table factor, so that conditions 
3 are also fulfilled; the con- 
ditions 4 are of themselves ful- 
filled. 

The equations for the 
stream function corresponding 
to equations (18) and (19) are 



!Pi = - 



-Jmi 



jjjnJL I b * bin (C + b ) 

t J! f 
bdln (£ + c) b (b — cl 



bjb-c) 



1 C[C + e) 



f (C + e) ^ C+e 
(6— c) In + (b — c)d In (f+ c ) 



? 



with d - - 1 - 2 

- a- 

The underlined terms are real, 
hence may "be omitted. 
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as before. ia identical with the follow- 
ing! 



bd 



u 



■6d 

•26 (6 



la(C_+ c)- 



in Hr+b 



IV. NUMERICAL EVALUATION 



The components of the ve- 
locity in the tunnel are: 



u 

U 



_d_ 
d x 



t(t + c) a*(£ + e) 



(6 — c)d 



j tog ln Tl 



6 — c 



6 — c 
"C + c* 



d a: 



"C/"*" D" 

'CI" 1 " # 

. . . (20) 



Only 



-(6- C )| 



In 



3L 
U 



and 



f- + • 



2. 
U 



at the profile 



In 



+ -> r ln6 



In 



edge are "being indicated. To 
calculate the differential quo- 

TT and IT 

ceed to the polar coordinates 
of the circular plane and 
that at the profile 



(fH} tl . nt . *u 



we pro- 



r ,cp 
note 



In 6 



i — ») (l-d){ ! 



In (f + 6) 



a 8 In 6 



In 



c (a 2 + c £ 

f).4 

; c J 



5 0 t 
d x 

\dyjr 



For 



The "boundary condition for 
\|/ = 0 at. the circle r » a was 
complied with "by deducting from 
each term the expression ob- 
tained if t " T 18 replaced 



[4*1 

ida:/r = a 

we have 



and 



- 0: 

r»a 

^ /r=a 
_ ( d(p ) r =a [ 
d y /r = a 



dxjr= a 
if.] 



since 



+ 6 



-c-H = sui y; hrH = 1 H cos 

\o#/r=a'. ■ -a \o.i//r=o a 

^-)a£COrdingto (19) has the form: 



by 



or 



t 1 



as reflected in 



I 

the first seven lines of the 
formula for V lt the last two 
*tabeing the remaining boundary 
conditions in order. If the 
procedure of the potential func- 
tion had been followed, it would 
have largely afforded the same 
terms as these, except with a dif- 
ferent prefix in some parts. 
The first term, for instance, 



U [dip Jr—a 



— (6 — c) d 



2i 



(b-' 



Vie 1 * 



(C + b) (C + e) 



2ie">" 



(c-6)(l-rf)| 2i ^ r 



(f + 6) (C + c) 



/ 2 t V y In (£ + 6) 

2 £ e^ln (£ + 6)_ 
" (t+c) 



(£+c) 2 JJr = a* 
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the maximum velocity at the profile in the tunnel is: 



U 
with 



Ufflax = 2 + 0 M s (21) 



e = l + f l + V) n + ^ I ~ ^ /i - *| sin_1 f 



-( l -£)-(*-£)J} 



Figure 3 shows the parameter 0 plotted against the 
displacement ratio. Por the case of the circle in free 

air a-5>0, 4 = 0, we get 0 = 2., as it already must be 
h 6 

the case according to Rayleigh's own calculated second 

approximation at the circle. For the other extreme 

a— »b, 4.= 0.5 (very long plate of thickness i n tun- 

h 2 

nel of height h) we get 0 = 3. In this case the veloc- 
ity "between profile and wall is constant. It follows 
from the requirement that the stream density is twice as 
great as at infinity: 



or 



u p = 2 U p 

U _ 3 Pot 
U p 

The expansion -222. according to the Mach number M = 



gives 



hence 0=3. 



H. = 2 + 3 M 3 + 
U 



The second curve plotted in figure 3 indicates the 
value of 0 in the narrowest area on the tunnel wall, 
when putting: 

33, = ^ik + 0 M s 

U U 
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where u ik is the velocity at the wall in incompressible 

flow. In figure 4 the entry of sonic velocity on the pro- 
file (critical Machmumber of flow) is plotted against the 
displacement ratio. '.'The-, critical Mach number' is obtained 

from equation (21) by putting - &1 ( a * = critical 

velocity ) . 

Figures 5 and 6 show the velocity distributions on 
the profile and in the narrowest part of the tunnel plot- 
ted against the Mach number of the flow. Theprofile 3 

of figure 1 a vith thickness ratio f = 0.761 and dis- 

t 

placement ratio ~ = 0.356, was chosen as cylinder. 



Translation by J. 
National Advisory 
for Aeronautics 



Vani er , 
Committ ee 
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Figure la.g Group of cylinders explored in tunnel, whose maximum velocity- 
is twice that of the approaching stream. 




Figure lb.- |' , n' plane. 
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Figure 2.- The displacement against thickness, 
d profile thickness 
t profile chord 
h tunnel height 
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Figure 3.- Effect of compressibility Figure 4.- Critical Each number. 

on the additional velocity 
at profile and tunnel wall. 
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Figure 5.- Velocity distribution Figure 6.- Velocity distribution in 

along profile of the minimum section for 

thickness ratio (1,0.761 (dJD.356) profile with thickness ratio d.0.761 

t h (d,0.356). t" 

plotted against polar angle. h" 
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